Vortex states in iron-based superconductors with collinear antiferromagnetic cores 
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Magnetism in the FeAs stoichiometric compounds and its interplay with superconductivity in 
vortex states are studied by self-consistently solving the BdG equations based on a two-orbital 
model with including the on-site interactions between electrons in the two orbitals. It is revealed 
that for the parent compound, magnetism is caused by the strong Hund's coupling, and the Fermi 
surface topology aids to select the spin-density-wave (SDW) pattern. The superconducting (SC) 
order parameter with s± = Ao cos^) cos(k y ) symmetry is found to be the most favorable pairing 
for both the electron- and hole-doped cases, while the local density-of-states (LDOS) exhibits the 
characteristic of nodal gap for the former and full gap for the latter. In the vortex state, the emer- 
gence of the field-induced SDW depends on the strength of the Hund's coupling and the Coulomb 
repulsions. The field-induced SDW gaps the finite energy contours on the electron and hole pocket 
sides, leading to the dual structures with one reflecting the SC pairing and the other being related 
to the SDW order. These features can be discernable in STM measurements for identifying the 
interplay between the field-induced SDW order and the SC order around the core region. 

PACS numbers: 74.20.Mn, 74.25.Ha, 74.25.Jb, 74.72.Bk 



I. INTRODUCTION 



The recently discovered iron arsenide superconduc- 
torSfii 2 -^^ which display superconducting transition 
temperature as high as more than 50K, appear to share a 
number of general features with high-T c cuprates, includ- 
ing the layered structure and proximity to a magnetically 
ordered stateii&i The accumulated evidences have sub- 
sequently established a fact that the parent compounds 
are generally poor metal and undergo structure and an- 
tiferromagnetic (AFM) spin-density-wave (SDW) transi- 
tions below certain temperatures Elastic neutron scat- 
tering experiments have shown the antiferromagnetic or- 
der is collinear and has a wavevector (n, 0) or (0, ir) in the 
unfolded Brillouin zone corresponding to a unit cell with 
only one Fe atom per unit cellp^ Either chemical dop- 
ing or/and pressure suppresses the AFM SDW instabil- 
ity and eventually results in the emergence of supercon- 
ductivity*^ The novel magnetism and superconducting 
properties in these compounds have been a great spur to 
recent researches ^^^li^^^i^ 2 ^ 

The relation between magnetism and superconductiv- 
ity and the origin of magnetic order have attracted signifi- 
cant attentions in the current research on FeAs supercon- 
ductors. Discrepancies exist in the experimental results, 
i.e., whether the superconductivity and antiferromag- 
netic order are well separated or they can coexist in the 
underdoped region of the phase diagram, and how they 
coexist if they happen to do so. For example, there is no 
overlap between those two phases in CeFeAsOi-^Fa^, 
while the coexistence of the two phases was observed 
in a narrow doping range in SmFeAsOi-^Fa— 2 -, and in 
a broader range in Bai_ a; K a ;Fe2As 2 11112 . Even for the 
same LaFeAsOi-^F^ system, different experiments dis- 
play conflicting results. It was reported that before the 



orthorhombic SDW phase is completely suppressed by 
doping, superconductivity has already appeared at low 
temperatures^, while it was also observed experimen- 
tally that superconductivity appears after the SDW is 
completely suppressed^. As for the origin of the SDW 
phase, two distinct types of theories have been proposed: 
local moment antiferromagnetic ground state for strong 
coupling^ and itinerant ground state for weak cou- 
pling] 15 i 16 ' 17 i 18 The detection of the local moment seems 
to question the weak coupling scenario, but the metallic- 
like (or bad metal) nature as opposed to a correlated 
insulator as in cuprates renders the strong coupling the- 
ories doubtableJ^ More recently, a compromised scheme 
was adopted: the SDW instability is assumed to result 
from the coupling of itinerant electrons with local mo- 
ment, namely, neither the Fermi surface nesting nor the 
local moment scenario alone is able to account for it. 20 

Although many research efforts have been already 
made to identify the existence of magnetic order and its 
origin as well as the relationship with superconductiv- 
ity, there have been fewer studies on vortex states of the 
systems. While the interplay between magnetism and su- 
perconductivity has been yet to be experimentally clari- 
fied, the superconducting critical temperature T c reaches 
its maximum value after the antiferromagnetic spin order 
is completely suppressed in the materials, indicating the 
competition nature between AFM SDW instability and 
superconductivity. At this stage, it is valuable and inter- 
esting to investigate vortex states in the family of FeAs 
compounds, mainly considering that the magnetic order 
may arise naturally when the superconducting order is 
destroyed by the magnetic vortex. Therefore, one can 
perform local tunneling spectroscopic probes in vortex 
states to understand profoundly the interplay between 
magnetic order and superconductivity. 
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In this paper, we investigate magnetism in the FeAs 
stoichiometric compounds, and the interplay between it 
and superconductivity upon doping in vortex states by 
self-consistently solving the BdG equations based on the 
two-orbital model with including the on-site interactions 
between electrons in the two orbitals. It is shown that 
for the parent compound, magnetism is caused by the 
strong Hund's coupling, and the Fermi surface topology 
aids to select the SDW ordering pattern. The SDW re- 
sults in the pseudogap-like feature at the Fermi level in 
the LDOS. It is found that the SC order parameter with 
s± = Aq cos(k x ) cos(ky) symmetry is the most favorable 
pairing at both the electron- and hole-doped sides, while 
the LDOS exhibits the characteristic of nodal gap for the 
former and full gap for the latter. In the vortex states, 
the emergence of the field-induced SDW order depends 
heavily on the strength of the Hund's coupling and the 
Coulomb repulsions. The coexistence of the field- induced 
SDW order and SC order around the core region is real- 
ized due to the fact that the two orders emerge at differ- 
ent energies. The corresponding LDOS at the core region 
displays a kind of dual structures, with one reflecting the 
SC pairing and the other being related to the SDW order. 

The paper is organized as follows. In Sec. II, we in- 
troduce the model Hamiltonian and carry out analytical 
calculations. In Sec. Ill, we present numerical calcu- 
lations and discuss the results. In Sec. IV, we make 
remarks and conclusion. 



II. THEORY AND METHOD 



stands for the vector potential in the Landau gauge and 
$o = hc/2e is the superconducting flux quantum. The 
hopping integrals are chosen as to capture the essence of 
the density function theory (DFT) results. 23 Taking the 
hopping integral between the d yz orbitals |<i| = 1 as the 
energy unit, we have, 
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(3) 

Here, x and y denote the unit vector along the x and y 
direction, respectively. 

The second term accounts for the superconducting 
pairing. Considering that a main purpose here is to ad- 
dress the interplay between the SC and magnetism in the 
vortex state for the FeAs-based superconductors, we take 
a phenomenological form for the pairing interaction, 
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with Vij as the strengths of effective attractions. 

The third term represents the interactions between 
electrons,— 
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We start with an effective two-orbital models that 
takes only the iron d xz and d yz orbitals into account. 
By assuming an effective attraction that causes the su- 
perconducting pairing and including the possible inter- 
actions between the two orbitals' electrons, one can con- 
struct an effective model to study the vortex physics of 
the iron-based superconductors in the mixed state: 

H = Hq + H pa i r + H int . (1) 

The first term is a tight-binding model 

-V E C L,<T C *,a, CT , ( 2 ) 



which describes the electron effective hoppings between 
sites i and j of the Fe ions on the square lattice, including 
the intra- (tij^ aa ) and inter-orbital (tij, a ,p, ot 7^ 13) hop- 
pings with the subscripts a, (3 (a, (/?) = 1,2 for xz and yz 
orbital, respectively) denoting the orbitals and a the spin. 
c\ aG creates an a orbital electron with spin a at the site 
i (i = (i x , i y )), and fx is the chemical potential. The mag- 
netic field is introduced through the Peierls phase factor 
e^w with Vij = ^ /;< A(r) • dr, where A = (-Hy, 0, 0) 
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which includes the intra-orbital (inter-orbital) Coulomb 
repulsion U (U ), the Hund's rule coupling J as well as 
the inter-orbital Cooper pairing hopping term J . 

After the Hartree-Fock decomposition of the on-site in- 
teraction term, one arrives at the Bogoliubov-de Genncs 
equations in the mean field approximation for this model 
Hamiltonian 
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where, 

-A 



(7) 



( u j,/3,<r)> w "a, CT ( u £/3,ir) are the Bogoliubov quasi- 
particle amplitudes on the j'-th site with corresponding 
eigenvalues E n . 

The pairing amplitude and electron densities are ob- 
tained through the following self-consistent equations, 25 

y 

^ij(i¥ : j)i ao > ~ "J / X U i,a,<J V i,a,g- v i,a,a u j,a,a) X 
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The electronic structure associated with the SDW 
and the vortex states, namely, the local density-of-states 
(LDOS), N(n,E) is calculated by 

N(r l ,E)= -Y J \K aA ?f'{E n -E) 



+K ad \ 2 f'(E n +E)}. 



(9) 



where, / (E) is the derivative of the Fermi-Dirac distri- 
bution function with respect to energy. 

In numerical calculations, the undoped case is de- 
termined by the equality of the area enclosed by the 
electron- and hole-pocket in the unfolded Brillouin zone, 
which leads to n h = J2i,a,a( n c<r) / ( N x N v) = 2 - Tlie 
Coulomb interactions U,U , J, J are expected to satisfy 
the conventional relation U = U — 2 J and J = Jj 17 i 26 
In the literatures, U ~ 0.2W - 0.5W and J - 0.09W are 
expected i 15 ' 17 Here, W is the energy bandwidth, which 
is 12.4|ii| in our case. This gives rise to J ~ \ti^ and 
U ~ 2.2 J — 5.5 J. We have found numerically that the 
results presented here are not subject to the qualitative 
changes in the intermediate coupling range U ~ 3 J ~ 4 J, 
where the ground state is an AFM meta L 28 i 29 In the fol- 
lowing, the typical result with U — 3.5J will be pre- 
sented. We take Vij — for the normal state. In the SC 
state, Vij is chosen to give a short coherence length of a 
few lattice spacing being consistent with experiments. 30 
We use Vij = V = 2.0, fi = 1.75 (fi = 1.13), which gives 



rise to the filling factor n = J^i a cr( n a,a)/(N x N y ) = 2.2 
(n = 1.8) and the coherent peak of the SC order pa- 
rameter in the DOS being at A max ~ 0.4. Thus, we 
estimate the coherence length £o ~ EFa/\A max \^- 1 ~ 4a. 
Due to this short coherence length, presumably the sys- 
tem will be a type-II superconductor. The unit cell with 
size N x x N y = 40 x 20 and the number of such unit cells 
M x x My = 10 x 20 are used in the numerical calculations. 
In view of these parameters, we estimate the upper criti- 
cal field B C 2 ~ 130T. Therefore, the model calculation is 
particularly suitable for the iron-based type-II supercon- 
ductors such as CaFei-^CozAsF, Euo.7Nao.3Fe2As2 and 
FeTei-zSz , where the typical coherence length £ de- 
duced from the experiments is of a few lattice spacing 30 
and the upper critical field achieves as high as dozens of 
Tesla. 



III. RESULTS AND DISCUSSION 

A. SDW phase in the absence of the magnetic field 

In the absence of a magnetic field and pairing term, we 
obtain the collinear AFM SDW at the half filling. Fig. 
1(a) is the typical result with J — 0.96 for the real space 
distribution of the moment M, defined as Mi = J2 a M ijQ , 
with Mi }Ct — \{{ni t<Xt -\) — (^i,a,j.)) being the spin order 
defined on the a orbital. As can be seen in Fig. 1(a), the 
real space distribution of Mj antiferromagnetically alines 
along the x direction but ferromagnetically along the y 
direction. In Fig. 1(b) Fourier transformation of Mi 
gives an SDW order with wave vector Q = (jr, 0), which 
is consistent with experimental results in the undoped 
systems^ 2 - [For another initial input parameters, the 
degenerate configuration of Mi with wave vector (0, tt) 
can be obtained.] We note that the emergence of mag- 
netic order is heavily dependent on the Hund's coupling 
strength J. For J = 0, the magnetic ordered phase 
does not exist even with very large U and U . There- 
fore, magnetism itself is generated by the strong Hund's 
coupling, whereas the Fermi surface topology aids to se- 
lect the ordering pattern^ 3 This is a reminiscent of the 
spin freezing phase found in a three-orbital model rele- 
vant to transition metal oxide SrRuOs^ and may be a 
common feature respecting the magnetic order origin in 
multiple orbital systems involving the Hund's coupling 
interaction. 

In Fig. 1(c), we plot the LDOS NfaE) in the SDW 
state at a site with positive M,-, i.e., the spin-up site 
labeled A in Fig. 1(a). The electronic structure in 
the SDW state displays a clear pseudogap-like feature 
with a heavily depressed but nonvanishing density-of- 
state (DOS) at the Fermi energy, pointing to the metal- 
lic magnetic ordered state. The magnetic order derived 
pseudogap-like feature is consistent with the experimen- 
tal observation of partial gaps in the SDW state of the 
parent compounds 34 and may account for the pseudogap 
feature in several experiments^ 5 -^ 
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The pseudogap feature comes from a fact that when 
the SDW order with the wave vector Q is involved, there 
will be gaps on those parts of the Fermi surface which 
are best connected by the wave vector Q, while those 
who are not connected by the wave vector Q remain un- 
touched, leading to the partial gaps in the SDW states 
of the parent compounds^ We make this point more 
clearly in Fig. 1(d), in which the spectral weight distri- 
bution I(k) = f^-Aw A(k, w)dui is shown. Here, A(k, u>) 
is the single-particle spectral function and Aw is an in- 
tegration window. As shown, both the electron and hole 
Fermi pockets are partially gaped. 

B. Configuration of the order parameters 

In the search for the most favorable pairing symme- 
try, we consider all possible singlet pairings, including 
the extended s- and ci-wave symmetries, between the 
nearest, next-nearest, and third-nearest neighbor (NN, 
NNN, TNN) sites, as shown in Figs. 2(a)-2(c). The 
pairing amplitude of the s-wave symmetry has the same 
sign along the x and y directions for the NN or TNN 
sites pairing, resulting in the k dependent pairing form 
A(fc) = Ao[cos(fc a ;)-|-cos(fcj / )] for the NN sites pairing and 
A(fc) = A [cos(2fc a; )+cos(2fc y )] for the TNN sites pairing, 
respectively; and the same sign along the x = y and x = 
— y directions for the NNN sites pairing, resulting in the 
k dependent pairing form A(fc) = Ao[cos(fc x ) cos(fc J/ )]. 
The d-wave pairing, on the other hand, has amplitude 
+Ao along the x direction and — Aq along the y direc- 
tion for the NN or TNN sites pairing, resulting in the k 
dependent pairing form A(fc) = A [cos(k x ) — cos(fc y )] for 
the NN sites pairing and A(fc) = Ao[cos(2fc x ) — cos(2fc a )] 
for the TNN sites pairing, respectively; and +Ao along 
x = y direction and — Ao along x = —y direction for the 
NNN sites pairing, resulting in the k dependent pairing 
form A(fc) = Ao[sm(k x ) sm(ky)]. 

The introduction of pairing interaction suppresses the 
SDW order completely on both the electron- and hole- 
doped sides, and leads to the homogeneous SC order in 
real space. We carry out extensive calculations, and find 
that in the reasonable doping range the most favorable 
pairing symmetry is the intra-orbital pairing between 
NNN sites, 

^i,i+s+y,QQ = Ai^i—x — y.aa — 

Ai^i+x — y^aa — Ai.i— x +y,aa — Ai_ aa: (10) 

which leads to the s-^-wave pairing s± = 
A cos(fc x ) cos(fc y ), being consistent with that ob- 
tained before i 17 ' 38 i 39 Then, the superconducting order 
parameter Aj is expressed as, 

Ai = i^A i>aa . (11) 

a 

For the choice of V = 2.0 in this paper, one gets the 
amplitude A, ~ 0.12 for the s± SC order . 



C. Vortex states 

When a magnetic field is applied, the SC order param- 
eter around the vortex core is suppressed, so that the 
system may be driven into a vortex state. We find that 
there exists a critical Hund's coupling value J c separat- 
ing the regimes of two kinds of vortex states associated 
respectively with and without the field-induced SDW or- 
der. In the following, we address these two regimes in 
detail. 



1. The vortex state without the field-induced SDW order 

The vortex state without the field induced SDW order 
is stable when J is less than J ce = 0.9 on the electron- 
doped side with n = 2.2 and less than J c h = 1.25 on 
the hole-doped side with n = 1.8, respectively. Typical 
results on the nature of the vortex state are displayed in 
Fig. 3 for n — 2.2 with the Hund's coupling J = 0.85, 
for which no magnetic order is induced. As shown in Fig. 
3(a), each unit cell accommodates two superconducting 
vortices each carrying a flux quantum hc/2e. The SC 
order parameter Aj vanishes at the vortex core center 
and recovers its bulk value at the core edge with radium 
£i on the scale of coherent length £q- 

In Figs. 3(b) and 3(c) we plot the LDOS as a function 
of energy at the vortex core center in the absence of the 
field-induced magnetic order for electron-doped case with 
n = 2.2 and hole-doped case with n = 1.8, respectively. 
For comparison, we have also displayed the LDOS at the 
midpoint between two nearest neighbor vortices along the 
x direction, which resembles that for the bulk system. As 
seen from Figs. 3(b) and 3(c), when J = 0.85 for which 
no local SDW order is induced, the LDOS at the core 
center shows a single resonant peak within the SC gap 
edge for both the electron- and hole-doped cases, which is 
similar to that reported by other authors for the cuprates 
high-T c superconductors in the vortex stated However, 
the differences are obvious with respect to the position of 
the resonant peak and the line shape of the bulk LDOS 
between the electron- and hole-doped cases in despite of 
the same SC pairing symmetry considered here. More 
specifically, for the electron-doped case, the position of 
the in-gap resonant peak is almost at the Fermi level and 
the bulk system exhibits the T^-shaped LDOS curve, the 
typical characteristics which indicate a nodal SC gap. 
However, for the hole-doped case, the resonant peak de- 
viates from the Fermi level to a higher energy and the 
bulk system exhibits the [/-shaped LDOS curve, from 
which the conclusion for a full SC gap can be made. 

The notable differences can be qualitatively under- 
stood as follows: The Fermi surface of the FeAs supercon- 
ductors consists of hole Fermi surfaces around the T-point 
at (k x , k y ) = (0, 0) forming the hole-pocket, and the elec- 
tron Fermi surfaces around the M\ t 2 points at (tt, 0) and 
(0, 7r) forming the electron-pocket, respectively. Both 
Fermi pockets change their size upon doping as depicted 
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in Fig. 3(d), where only the relevant electron-pockets 
are displayed. The size of the electron-pocket enlarges 
and approaches to the nodal line of the s± SC gap with 
electron doping while shrinks and deviates from the nodal 
line with hole doping. Thus the low energy quasiparticles 
in the SC phase show the nodal behavior in the electron- 
doped system and nodclcss behavior in the hole-doped 
system. This may explain the discrepancy observed 
in experiments concerning the pairing symmetry, where 
the conclusion for the nodal gap were obtained on the 
electron-doped LnFeAsO (Ln stands for the rare earth el- 
ements) sample a 41 ' 42 ! 43 and a dominant full gap feature 
was found on the hole-doped (Ba,Sr)i_. c K, E Fe2As2 sys- 
tem o 34 ^ 44 in the measurement of the thermal and trans- 
port properties. 



2. The vortex state with the field-induced SDW order 

As J increases to about J ce = 0.9 on the electron- 
doped side with n = 2.2 and J c h = 1-25 on the hole- 
doped side with n = 1.8, the SDW order is induced 
around the core region. Fig. 4 displays the vortex struc- 
ture with J — 0.96 for the electron-doped case, where 
the local magnetic order is induced around the vortex 
core, as shown in Fig. 4(b) which presents the spatial 
distribution of the local SDW order as defined in Sec. 
Ill A. As seen in Fig. 4(a), the vortex core expands fur- 
ther with a radium £2 compared with that in Fig. 3(a). 
Meanwhile, the maximum strength of Mj appears at the 
vortex core center and decays with a scale of £2 to zero 
into the superconducting region, depicting a competition 
nature between SC and magnetic orders as observed in 
experiment 4 ^. 

In this case, it is shown that there is no obvious split- 
ting of the in-gap bound state peak though the peak in- 
tensity is suppressed heavily in the LDOS for both the 
electron- and hole-doped cases, as shown in Figs. 4(c) 
and 4(d). In addition to the in-gap state peak, an ad- 
ditional peak structure below the Fermi energy appears 
for the electron-doped case, while two peaks situate sep- 
arately around and above the Fermi energy for the hole- 
doped case. These features are dramatically different 
from the high-T c cuprates in vortex states with the field- 
induced antiferromagnetic order, where the in-gap reso- 
nant peak of the core bound state is split into two peaks 
sitting symmetrically about the Fermi energy4£ 

Simply, one can analyze the present vortex state in the 
following way: there are two factors that play a role in the 
physics around the vortex core region. One is the pure SC 
of vortex state without SDW in the magnetic field, while 
the other is the SDW state without the SC order for the 
doped case. It is known that doping destroys the nest- 
ing property between parts of the Fermi surface on the 
electron- and hole-pocket, due to the size change of the 
hole- and electron-pocket upon doping compared with 
that in the undoped case. However, as depicted in Figs. 
5(c) and 5(d), the SDW wave vector Q now connects the 



finite energy contour for the doped case, resulting in the 
gap-like feature below (above) Fermi energy in the LDOS 
shown in Fig. 5(a) (Fig. 5(b)) for electron (hole) doped 
case. Combination of the in-gap resonant peak in the 
vortex state without SDW and the SDW-induced finite 
energy gap feature produce a kind of dual structures of 
the LDOS for the finite doped case, ie., the in-gap bound 
state peak reflecting the SC pairing and the other peak 
structure being related to the SDW order. 

IV. REMARKS AND CONCLUSION 

Clarification of the interplay between magnetism and 
superconductivity is a key step toward the understanding 
of the underlying physics of the Fe-based high-T c super- 
conductors. Although the competition nature between 
them has been identified in both classes of materials, 
some still show a coexistence of ther m 11 ! 12 ] 22 ! 45 Com- 
petition between the AFM SDW and SC is natural in 
FeAs compounds when one considers that both originate 
from the multiple Fe d conduction bands, but quest for 
the mechanism of the coexistence of them is shown to 
be more challenged. Recently, an incommensurate SDW 
state with wave vector Q = Q ± q has been proposed 
to account for the coexistence of the AFM SDW and su- 
perconductivity at finite dopin g 47 i 48 . In such a state, the 
mismatch between the electron and hole Fermi pocket 
at finite doping is compensated by the incommensurate 
wave vector q, leading to the inferior "nesting" between 
the electron and hole Fermi pocket and allowing for the 
coexistence of magnetism and superconductivity. This 
mechanism only works at the doping level near the AFM 
instabilities, where the mismatch between the electron 
and hole Fermi pocket is small. Here, we show that the 
field-induced SDW at the doping level being far from the 
AFM instabilities is commensurate with the same wave 
vector Q as in the undoped case, and it gaps the finite 
energy contour on the electron and hole pocket sides. 
Therefore, at optimal doping, the field-induced SDW and 
the SC order around the core region may associate with 
the DOS at different energies, allowing them to coexist. 

In conclusion, we have studied magnetism in the FeAs 
stoichiometric compounds and the interplay between it 
and superconductivity upon doping in the vortex state 
by self-consistently solving the BdG equations based on 
the two-orbital model including the on-site interactions 
between electrons in the two orbitals. It has been shown 
that for the parent compound, magnetism is caused 
by the strong Hund's coupling, and the Fermi surface 
topology aids to select the SDW ordering pattern. The 
SDW results in the pseudogap-like feature at the Fermi 
level in the LDOS. The SDW order is completely sup- 
pressed upon the introduction of the SC interaction. 
We have also found that the SC order parameter with 
s± = A cos(k x ) cos(ky) symmetry is the most favorable 
pairing at both the electron- and hole-doped sides, while 
the LDOS exhibits the characteristic of nodal gap for 
the former and full gap for the latter. In vortex states, 



6 



the emergence of the field-induced SDW order depends 
heavily on the strength of the Hund's coupling and the 
Coulomb repulsions, while the coexistence of the field- 
induced SDW order and SC order around the core region 
is realized due to the fact that the two orders emerge at 
different energies. The LDOS at the core region for the 
vortex state with SDW displays the dual structures with 
one reflecting the SC pairing and the other being related 
to the SDW order. These features can be discernable 
in the STM measurements for identifying the interplay 
between the field-induced SDW order and the SC order 
around the core region. 

V. ACKNOWLEDGEMENT 

This project was supported by National Natural Sci- 
ence Foundation of China (Grant No. 10525415 and No. 
10904062), the Ministry of Science and Technology of 
Science (Grants Nos. 2006CB601002, 2006CB921800), 
a GRF grant of Hong Kong (HKU7055/09P), the 
China Postdoctoral Science Foundation (Grant No. 
20080441039), and the Jiangsu Planned Projects for 
Postdoctoral Research Funds (Grant No. 0801008C). 

Note added- Alter completing this work, we get to know 
a related work done by X. Hu, C. S. Ting, and J. X. 
Zhuia. 




FIG. 1: (Color online) (a) The real space distribution of the 
moment Mi. (b) The Fourier transformation of Mi. (c) The 
LDOS at the site labeled as A in Fig. 1(a) in the SDW state, 
and (d) spectral weight distribution in the SDW state (see 
text) . 



Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 
Am. Chem. Soc. 130, 3269 (2008). 

X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and D. 

F. Fang, Nature 453, 761 (2008). 

Z.-A. Ren, G.-C. Che, X.-L. Dong, J. Yang, W. Lu, W. Yi, 
X.-L. Shen, Z.-C. Li, L.-L. Sun, F. Zhou, and Z.-X. Zhao, 
Europhys. Lett. 83, 17002 (2008). 

G. F. Chen, Z. Li, D. Wu, G. Li, W. Z. Hu, J. Dong, P. 
Zheng, J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 100, 
247002 (2008). 

C. Wang, L. Li, S. Chi, Z. Zhu, Z. Ren, Y. Li, Y. Wang, X. 
Lin, Y. Luo, S. Jiang, X. Xu, G. Cao, and Z. Xu, Europhys. 
Lett. 83, 67006 (2008). 

C. de la Cruz, Q. Huang, J. W. Lynn, J. Li, W. Ratcliff II, 
J. L. Zarestky, H. A. Mook, G. F. Chen, J. L. Luo, N. L. 
Wang, and P. Dai, Nature 453, 899 (2008). 
J. Dong, H. J. Zhang, G. Xu, Z. Li, G. Li, W. Z. Hu, D. 
Wu, G. F. Chen, X. Dai, J. L. Luo, Z. Fang, and N. L. 
Wang, Europhys. Lett. 83, 27006 (2008). 

H. -H. Klauss, H. Luetkens, R. Klingeler, C. Hess, F. J. 
Litterst, M. Kraken, M. M. Korshunov, I. Eremin, S.-L. 
Drechsler, R. Khasanov, A. Amato, J. Hamann-Borrero, 
N. Leps, A. Kondrat, G. Behr, J. Werner, and B. Buchner 
Phys. Rev. Lett. 101, 077005 (2008). 

H. Takahashi, K. Igawa, K. Arii, Y. Kamihara, M. Hirano, 
and H. Hosono, Nature (London) 453, 376 (2008). 
J. Zhao, Q. Huang, C. de la Cruz, S. Li, J. W. Lynn, Y. 
Chen, M. A. Green, G. F. Chen, G. Li, Z. Li, J. L. Luo, 
N. L. Wang, P. Dai, Nature Materials 7, 953 (2008). 



g c o 



(-). 



+ (") 



-O D 



V.;) 



+(->. 



(a) 



(b) 



(c) 



FIG. 2: (Color online) The SC order parameter configuration 
in the real space between the NN sites pairing (a), the NNN 
sites pairing (b), and the TNN sites pairing (c), respectively. 
The minus signs in the parentheses correspond to the d-wave 
pairings, otherwise the s-wave pairings. 



H. Chen, Y. Ren, Y. Qiu, Wei Bao, R. H. Liu, G. Wu, T. 
Wu, Y. L. Xie, X. F. Wang, Q. Huang, and X. H. Chen, 
Europhys. Lett. 85, 17006 (2009). 

T. Goko, A. A. Aczel, E. Baggio-Saitovitch, S. L. Bud'ko, 
P. C. Canfield, J. P. Carlo, G. F. Chen, P. Dai, A. C. 
Hamann, W. Z. Hu, H. Kageyama, G. M. Luke, J. L. Luo, 
B. Nachumi, N. Ni, D. Reznik, D. R. Sanchez-Candela, 
A. T. Savici, K. J. Sikes, N. L. Wang, C. R. Wiebe, T. 



7 




FIG. 3: (Color online) (a) The real space distribution of the 
SC order amplitude |Aj| without the field-induced SDW. The 
LDOS curves at the core center (black line) and for the bulk 
system (green line) with electron doping n = 2.2 (b), and hole 
doping n — 1.8 (c), respectively, (d) The electron-pocket in 
the iron-pnictides in the unfolded Brillouin zone — n < k x < 
tv, — n < k y < 7r. The solid and dashed curves correspond to 
electron doping with n = 2.2 and hole doping with n = 1.8, 
respectively. The dotted (red) lines mark the nodal lines at 
(7r/2, fey) and (k x ,n/2) for the s± = Ao cos^) cos(fe H ) order 
parameter. 



J. Williams, T. Yamamoto, W. Yu, and Y. J. Uemura, 
larXiv:0808.T425l (unpublished) . 

H. Luetkens, H.-H. Klauss, M. Kraken, F. J. Litterst, T. 
Dellmann, R. Klingeler, C. Hess, R. Khasanov, A. Am- 
ato, C. Baines, M. Kosmala, O. J. Schumann, M. Braden, 
J. Hamann-Borrero, N. Leps, A. Kondrat, G. Behr, J. 
Werner, and B. Buechner, Nature Materials 8, 305 (2009). 
Q. Si, and E. Abrahams, Phys. Rev. Lett. 101, 076401 
(2008). Z. Y. Weng, Physica E 41, 1281 (2009); C. Xu, 
M. Miiller, and S. Sachdev, Phys. Rev. B 78, 020501(R) 
(2008); C. Fang, H. Yao, W.-F. Tsai, J. Hu, and S. A. 
Kivelson, Phys. Rev. B 77, 224509 (2008); G.-M. Zhang, 
Y.-H. Su, Z.-Y. L u, Z.-Y. Weng, D.-H. Lee, and T. Xiang, 
larXiv:0809.3874l (unpublished) 

C. Cao, P. J. Hirschfeld, and H.-P. Cheng, Phys. Rev. B 
77, 220506(R) (2008); S. Yang, W.-L. You, S.-J. Gu, and 
H.-Q. Lin, Chin. Phys. B 18, 2545 (2009). 

D. J. Singh and M. H. Du, Phys. Rev. Lett. 100, 237003 
(2008); K. Kuroki, S. Onari, R. Arita, H. Usui, Y. Tanaka, 
H. Kontani, and H. Aoki, Phys. Rev. Lett. 101, 087004 
(2008); Q. Han, Y. Chen, and Z. D. Wang, Europhys. Lett. 
82, 37007 (2008); Q. Han and Z. D. Wang, New J. Phys. 
11, 025022 (2009); V. Cvetkovic, and Z. Tesanovic, Euro- 
phys. Lett. 85, 37002 (2009). 

Z.-J. Yao, J.-X. Li, and Z. D. Wang, New J. Phys. 11, 
025009 (2009); S. L. Yu, J. Kang, and J. X. Li, Phys. Rev. 
B 79, 064517 (2009). 




FIG. 4: (Color online) (a) The real space distribution of the 
SC order amplitude |Aj| in the presence of the field-induced 
SDW. (b) The real space distribution of the field-induced mo- 
ment Mi. The LDOS curves at the core center (black line) and 
for the bulk system (green line) with electron doping n — 2.2 
(c), and hole doping n = 1.8 (d), respectively. 



18 S. Raghu, X.-L. Qi, C.-X. Liu, D. J. Scalapino, and S.-C. 
Zhang, Phys. Rev. B 77, 220503(R) (2008). 

19 W. Z. Hu, G. Li, J. Dong, Z. Li, P. Zheng, G. F. Chen, 
J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 101, 257005 
(2008). 

20 J. Dai, Q. Si, J.-X. Zhu, and E. Abrahams, 
larXiv:0808.0305l (unpublishe d); S.-P. Kou, T. Li, and 
Z.-Y. Weng, [arXTv:081 1.41 11 1 (unpublished) . 

21 Y. Wan, Q.-H. Wang, EPL 85, 57007 (2009). 

22 R. H. Liu, G. Wu, T. Wu, D. F. Fang, H. Chen, S. Y. Li, 
K. Liu, Y. L. Xie, X. F. Wang, R. L. Yang, L. Ding, C. 
He, D. L. Feng, and X. H. Chen, Phys. Rev. Lett. 101, 
087001 (2008); A. J. Drew, F. L. Pratt, T. Lancaster, S. 
J. Blundell, P. J. Baker, R. H. Liu, G. Wu, X. H. Chen, 
I. Watanabe, V. K. Malik, A. Dubroka, K. W. Kim, M. 
Rossle, and C. Bernhard, Phys. Rev. Lett. 101, 097010 
(2008). 

23 G. Xu, W. Ming, Y. Yao, X. Dai, S.-C. Zhang and Z. Fang, 
Europhys. Lett. 82, 67002 (2008); Z. P. Yin, S. Lebegue, 
M. J. Han, B. P. Neal, S. Y. Savrasov, and W. E. Pickett, 
Phys. Rev. Lett. 101, 047001 (2008). 

24 A. M. Oles, G. Khaliullin, P. Horsch, and L. F. Feiner, 
Phys. Rev. B 72, 214431 (2005). 

25 Y. Chen, Z. D. Wang, J.-X. Zhu, and C. S. Ting, Phys. 
Rev. Lett. 89, 217001 (2002). 

26 P. Werner, E. Gull, M. Troyer, and A. J. Millis, Phys. Rev. 
Lett. 101, 166405 (2008). 

27 K. Haule, and G. Kotliar, New J. Phys. 11, 025021 (2009). 

28 J. Lorenzana, G. Seibold, C. Ortix, and M. Grilli, Phys. 
Rev. Lett. 101, 186402 (2008). 

29 R. Yu, K. T. Trinh, A. Moreo, M. Dagho fer, J. A. Riera, 
S. Haas, and E. Dagotto. [arXiv:0812. 28941 (unpublished - ) . 



8 





2 -1 



Bit, I 

(d) 


1 ; 


J 




V 

Q - 




o 


(C 




fib 





-0.5 0.0 0.5 1.0 



FIG. 5: (Color online) LDOS curves in a bulk system simu- 
lating the SDW state around the core region for n = 2.2 (a), 
and n — 1.8 (b), respectively. The wave vector Q of the field- 
induced SDW connects the finite energy contour (red line) for 
7i = 2.2 (c), and n = 1.8, respectively (Black curves denote 
the Fermi surface at corresponding doping level). 



30 S. Takeshita, R. Kadono, M. Hiraishi, M. Miyazaki, A. 
Koda, S. Matsuishi, and H. Hosono, Phys. Rev. Lett. 103, 
027002 (2009); Y. Qi, Z. Gao, L. Wang, D. Wang, X. 
Zhang, and Y. Ma, New J. Phys. 10, 123003 (2008); Y. 
Mizuguchi, F. Tomioka, S. Tsuda, T. Yamaguchi, and Y. 
Takano, Appl. Phys. Lett. 94, 012503 (2009). 

31 Y. D. Zhu, F. C. Zhang, and M. Sigrist, Phys. Rev. B 51, 
1105 (1995). 

32 Y. Chen, J. W. Lynn, J. Li, G. Li, G. F. Chen, J. L. Luo, 
N. L. Wang, P. Dai, C. de la Cruz, and H. A. Mook, Phys. 
Rev. B 78, 064515 (2008). 

33 I. I. Mazina and J. Schmalian, Physica C 469, 614 (2009). 

34 W. Z. Hu, Q. M. Zhang, N. L. Wang, Physica C 469, 545 
(2009). 

35 H.-Y. Liu, X.-W. Jia, W.-T. Zhang, L. Zhao, J.-Q. Meng, 
G.-D. Liu, X.-L. Dong, G. Wu, R.-H. Liu, X.-H. Chen, 
Z.-A. Ren, W. Yi, G.-C. Che, G.-F. Chen, N.-L. Wang, 

G. -L. Wang, Y. Zhou, Y. Zhu, X.-Y. Wang, Z.-X. Zhao, 
Z.-Y. Xu, C.-T. Chen, and X.-J. Zhou, Chin. Phys. Lett. 
25, 3761 (2008); L. Zhao, H.-Y. Liu, W.-T. Zhang, J.-Q. 
Meng, X.-W. Jia, G.-D. Liu, X.-L. Dong, G.-F. Chen, J.-L. 
Luo, N.-L. Wang, W. Lu, G.-L. Wang, Y. Zhou, Y. Zhu, 
X.-Y. Wang, Z.-Y. Xu, C.-T. Chen, and X.-J. Zhou, Chin. 
Phys. Lett. 25, 4402 (2008). 

36 M. C. Boyer, K. Chatterjee, W. D. Wise, G. F. Ch en, J. L. 
Luo, N. L. Wang, E. W. Hudson. larXiv:0806.4400l (unpub- 
lished); T. Sato, S. Souma, K. Nakayama, K. Terashima, 
K. Sugawara, T. Takahashi, Y. Kamihara, M. Hirano, and 

H. Hosono, J. Phys. Soc. Jpn. 77, 063708 (2008); Y. Ishida, 
T. Shimojima, K. Ishizaka, T. Kiss, M. Okawa, T. To- 



gashi, S. Watanabe, X. Y. Wang, C. T. Chen, Y. Kami- 
hara, M. Hirano, H. Hosono, and S. Shin, Phys. Rev. B 
79, 060503(R) (2009). 

37 D. Hsieh, Y. Xia, L. Wray, D. Qian, K.K. Gomes, A. Yaz- 
dani, G. F. Chen, J. L. Luo, N. L. Wang, and M. Z. Hasan, 
larXiv:0812.2289l (unpublished) . 

38 I. I. Mazin, D. J. Singh, M. D. Johannes, and M. H. 
Du, Phys. Rev. Lett. 101, 057003 (2008); K. Seo, B. A. 
Bernevig, and J. P. Hu, Phys. Rev. Lett. 101, 206404 

(2008) . 

39 F. Wang, H. Zhai, Y. Ran, A. Vishwanath, and D. H. 
Lee, Phys. Rev. Lett. 102, 047005 (2009); W.-Q. Chen, 
K.-Y. Yang, Y. Zhou, and F.-C. Zhang, Phys. Rev. Lett. 
102, 047006 (2009); A. V. Chubukov, Physica C 469, 640 

(2009) . 

40 Y. Wang, and A. H. MacDonald, Phys. Rev. B 52, R3876 
(1995); M. Franz, and Z. Tesanovic, Phys. Rev. Lett. 80, 
4763 (1998). 

41 G. Mu, X. Y. Zhu, L. Fang, L. Shan, C. Ren, and H. 
H. Wen, Chin. Phys. Lett. 25, 2221 (2008); L. Shan, Y. L. 
Wang, X. Y. Zhu, G. Mu, L. Fang, C. Ren, and H. H. Wen, 
Europhys. Lett. 83, 57004 (2008); C. Ren, Z. S. Wang, H. 
Yang , X. Y. Zhu L. Fang, G. Mu, L. Shan, and H. H. Wen, 
larXiv:0804.T726l (unpublished). 

42 H. -J. Grafe, D. Paar, G. Lang, N. J. Curro, G. Behr, J. 
Werner, J. Hamann-Borrero, C. Hess, N. Leps, R. Klin- 
geler, and B. Buechner, Phys. Rev. Lett. 101, 047003 
(2008); H. Luetkens, H.-H. Klauss, R. Khasanov, A. Am- 
ato, R. Klingeler, I. Hellmann, N. Leps, A. Kondrat, 
C.Hess, A. Kohler, G. Behr, J. Werner, and B. Buchner, 
Phys. Rev. Lett. 101, 097009 (2008); J. P. Carlo, Y. J. 
Uemura, T. Goko, G. J. MacDougall, J. A. Rodriguez, W. 
Yu, G. M. Luke, P. C. Dai, N. Shannon, S. Miyasaka, S. 
Suzuki, S. Tajima, G. F. Chen, W. Z. Hu, J. L. Luo, and 
N. L. Wang, Phys. Rev. Lett. 102, 087001 (2009). 

43 K. Matano, Z. A. Ren, X. L. Dong, L. L. Sun, Z. X. Zhao, 
and G. Q. Zheng, Europhys. Lett. 83, 57001 (2008); K. 
Ahilan, F. L. Ning, T. Imai, A. S. Sefat, R. Jin, M. A. 
McGuire, B. C. Sales, and D. Mandrus, Phys. Rev. B 78, 
100501 (R) (2008); F. Massee, Y. Huang, R. Huisman, S. de 
Jong, M. S. Golden, and J. B. Goedkoop, rarXiv:0812.4539l 
(unpublished) . 

44 G. Li, W. Z. Hu, J. Dong, Z. Li, P. Zheng, G. F. Chen, 
J. L. Luo, and N. L. Wang, Phys. Rev. Lett. 101, 107004 

(2008) ; G. Mu, H. Luo, Z. Wang, L. Shan, C. Ren, and 
Hai-Hu Wen, Phys. Rev. B 79, 174501 (2009). 

45 D. K. Pratt, W. Tian, A. Kreyssig, J. L. Zarestky, S. 
Nandi, N. Ni, S. L. Bud'ko, P. C. Canfield, A. I. Gold- 
man, and R. J. McQueeney, Phys. Rev. Lett. 103, 087001 

(2009) . 

46 I. Maggio-Aprile, Ch. Renner, A. Erb, E. Walker, and 0. 
Fischer, Phys. Rev. Lett. 75, 2754 (1995); J.-X. Zhu, and 
C. S. Ting, Phys. Rev. Lett. 87, 147002 (2001). 

47 V. Cvetkovic and Z. Tesanovic, Phys. Rev. B 80, 024512 
(2009). 

48 A. B. Vorontsov, M. G. Vavilov, and A. V. Chubukov, 
Phys. Rev. B 79, 060508(R) (2009). 

49 X. Hu, C. S. Ting, and J. X. Zhu, Phys. Rev. B 80, 014523 
(2009). 



